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Abstract. We present GMRT observations of the Galactic microquasar GRS 1915+105 at 1.28 GHz for 8 days 
from 2001 June 18 to July 1. We have seen several isolated radio flares of varying magnitudes (20 - 50 mjy) 
and durations (6 - 35 min) and we model them as due to adiabatically expanding synchrotron emitting clouds 
(synchrotron bubbles) ejected from the accretion disk. By applying this model, we provide a new method to 
estimate the electron power-law index p, hence the spectral index, from single frequency radio observations. 
This method does not require correction for the optical time delay effects which may be important in the case 
of optically thick radio emission. Using our estimated value of p and simultaneous multiwavelength data from 
literature, we have calculated the time of ejection of the synchrotron plasma and the time delays at different 
observed frequencies. Our estimates are in good agreement with the observed time delays. 
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1. Introduction 

The Galactic X-ray transient source GRS 1915+105 was 
discovered in 1992 by the GRANAT satellite (Castro-Tirado 
et al. 1992). It has now been established as a black hole 
binary with the black hole mass of 14±4Mq and the com- 
panion K-M III star of mass 1.2±O.2M0 (Greiner et al. 
2001). The source was shown to possess relativistic out- 
flows during its radio outbursts (Mirabel & Rodriguez 
1994). Extensive monitoring of GRS 1915+105 in the ra- 
dio and X-ray band showed that the overall radio emission 
is correlated with its X-ray properties (eg: Harmon et al. 

1997) . The radio emission can be broadly classified into 
three classes; (i) the relativistic superluminal radio jets of 
flux density ~ 1 Jy with decay time-scales of several days 
(Fender et al. 1999), (ii) the baby jets of 20 - 40 min du- 
rations with flux density of 20 — 200 mJy both in infrared 
(IR) and radio (Pooley & Fender 1997; Eikenberry et al. 

1998) and (iii) the plateau state with persistent radio emis- 
sion of 20 — 100 mJy for extended durations (Muno et al. 
2001). In the case of superluminal jets, the radio emission 
has steep spectra and are observed at large distances (400 
- 5000 AU) from the accretion disk (Fender et al. 1999; 
Dhawan et al. 2000). The radio emission at this distance 
is believed to be decoupled from the accretion disk (Yadav 
2001; Muno et al. 2001). In contrast, the other two classes 
of radio emission has flat spectra and occur close to the 
accretion disk (within a few tens of AU). Even though su- 



perluminal jets and baby jets are differentiated by their 
spectra, decay time scales, and the distance (from the ac- 
cretion disk) where the emission occurs, there is evidence 
for ejection of significant amount of relativistic material 
even during the baby jets (Fender et al. 1999). The ra- 
dio emission is believed to be due to synchrotron emission 
from relativistic electrons and the dominant decay mecha- 
nism is the adiabatic expansion losses (Mirabel et al. 1998, 
hereinafter M98). Feroci et al. (1999) have studied an iso- 
lated radio flare (~ 60 mJy at 15 GHz) in X-ray and ra- 
dio and present strong evidence for mass ejection (QPO 
disappears during the radio outburst) . Coordinated multi- 
wavelength studies have indicated strong disk-jet connec- 
tion (M98; Eikenberry et al. 1998; Yadav 2001). 

The time of the ejection of the radio emitting syn- 
chrotron plasma from the accretion disk is an impor- 
tant parameter to understand the disk— jet connection. 
Another important parameter is the power-law index (p) 
of the electron energy distribution in the plasma which 
is related to the delay time. Even though simultaneous 
multi- wavelength radio observations provide p, in the case 
of optically thick radio emission it is affected by the opti- 
cal depth time delay effects and may introduce large error 
in the estimate of time of ejection, particularly in the case 
of periodic emission like IR and radio baby jets. 

In this Letter we provide a new method, by applying 
the adiabatic expansion model to baby jets to estimate 
p from single frequency radio observations. This method 
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the optical depth effects. We apply this method to isolated 
baby jets seen in the radio observations of GRS 1915+105 
with the Giant Meterwave Radio Telescope (GMRT) at 
1.28 GHz and estimate p and the time of ejection of the 
synchrotron bubble. Using this p and the multiwavelength 
data from literature we also calculate the time delays at 
different observed frequencies and the time at which the 
synchrotron bubble is ejected. In Section 2 we present the 
observational details and basic results. The adiabatic ex- 
pansion model is given in Section 3. Discussions and con- 
clusions are presented in Section 4. 

2. Observations and Results 

The radio observations were carried out at 1.28 GHz with 
a bandwidth of 16 MHz using the GMRT (Swarup ct al. 
1991) on 2001 June 18, 22, 23, 27, 28, 29, 30 and July 
1. The flux density scale is set by observing the primary 
calibrator 3C286 or 3C48. A phase calibrator was ob- 
served before and after a 45 min scan on GRS 1915+105. 
The integration time was 32 s. The data recorded from 
GMRT has been converted to FITS and was analysed us- 
ing Astronomical Image Processing System (MPS). The 
lightcurve produced at 5 min interval is shown in Figure 
1 (top panel). GRS 1915+105 exhibited significant radio 
emission on all days, except on June 27. On June 18, the 
source was quite steady with the flux density ~ 10 mJy. 
However, on June 22 the source became bright (~ 50 mJy) 
and showed rapid drop in the flux density by a factor 
of two and slowly recovered to about the same level. On 
June 23, the flux density was steady at about 30 mJy. 
The source became very weak in radio (~ 5 mJy) on June 

27. However on 28, the flux density was ~ 5 mJy at about 
UTC 17 hour (beginning of observation) and reached grad- 
ually to 50 mJy at UTC 24 hour. When the observations 
resumed on 29, the source was "caught" at 70 mJy at UTC 
16 hour, and the flux started decaying slowly to the value 
of 10 mJy at UTC 24 hour. The source was relatively faint 
on June 30 and July 1. 

We have checked the radio data at shorter binning time 
of 1 min, which showed clear short term variations. Minor 
flares of amplitude upto 50 mJy have been seen on June 

28, 29, 30 and July 1, which are typical of baby jets (Figure 
3 and Table 1). The X-ray hardness ratio (5-12 keV/1.5- 
5 keV) from RXTE/ASM shown in the bottom panel of 
Figure 1 suggests that overall GRS 1915+105 remained in 
the low hard state. 

3. Adiabatic Expansion Model 

We propose that the minor flares (baby jets) seen in 
the GMRT observations are due to adiabatically expand- 
ing synchrotron emitting bubbles. Hjellming & Johnston 
(1988) present models for explaining the persistent radio 
emission from X-ray binaries by assuming adiabatic ex- 
pansion dominated by lateral motions. For flaring events, 
they consider spherical bubbles of relativistic plasma ex- 
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Fig. 1. Top panel: Radio light curve from GMRT at 1.28 
GHz, binned at 5 min interval. The observations on 2001 
June 18, 22, 23 and 27 is for ~ 30— 60 min, while rest of 
the days are of 8-10 hours. Bottom panel: RXTE/ASM 
hardness ratio (5-12 keV/1.5-5 keV) for the same duration. 

batically expanding sphere model (three dimensional ex- 
pansion) for the baby jets seen in the GMRT observations 
to understand their overall radio behavior. We assume 
that the expansion is uniform and the transition time from 
optically thick to thin state is negligible. 

Consider a spherical cloud of relativistic electrons with 
the energy distribution N{E)dE oc E~PdE, emitting syn- 
chrotron radiation. If this cloud is expanding adiabatically, 
the observed flux density as a function of time will be (van 
der Laan, 1966; Hjellming & Johnston, 1988) 



under optically thick conditions, and 
Se/Sp = (re/rp)-2p 



(1) 



(2) 

under optically thin conditions. Here Ss, Sp and are the 

observed flux densities at the start, peak and the end of 
the flare and rg , rp and are the corresponding radii of 
the bubble, p is the power-law index of electron energy 
distribution. We define the time sequence of synchrotron 
bubble as in Figure 2. is the time taken for the bubble 
to attain the flux density Sg, tp is the time taken since 
tg for the flux density to reach its peak value Sp (expan- 
sion under optically thick conditions) and is the time 
taken since tp by the bubble to decay to the flux level of 
Sc (expansion under optically thin conditions). Assuming 
uniform expansion, the profiles will be 



Sg[l + itp/ts)f 



(3) 



for the rising part of the light curve and 

S, = Sp{l + [U/itp+U)]}-^P (4) 

for the decay part. An interesting case is when the syn- 
chrotron bubble is isolated, where the rise, peak and decay 
of a single bubble is observed. The adiabatic equations (1) 
and (2) can be solved to obtain p 
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Fig. 2. Illustration of the time markers during the adia- 
batic expansion. The point on the left marks the ejection 
time. 



Thus by measuring the rise and decay time from isolated 
flares, the electron power-law index p can be obtained. For 
estimating p, we have considered only those flares where 
the rise and decay of single flare is clearly visible well 
above the noise level or steady emission. We fit the above 
equations (using the least square method) to such isolated 
flares seen in our data on 2001 June 28, 29, 30 and July 1. 
The best fit model for two such flares are given in Figure 

3. The estimated parameters are given in Table 1. Column 
1 gives the date of observation, Column 2 gives the esti- 
mated time (UT) of birth of the synchrotron bubble (see 
Section 4 for details). The duration and peak flux density 
of the flares are given in Columns 3 and 4. The electron 
power-law index (p) and the spectral index {a = (p— 1)/2) 
are given in Columns 5 and 6. 

4. Discussion and Conclusion 

We have presented the GMRT observations of the micro- 
quasar GRS 1915-1-105 at 1.28 GHz. Several isolated mi- 
nor flares which are typical of baby jets have been found. 
We have suggested that these flares are due to adiabati- 
cally expanding bubbles of synchrotron emitting plasma. 
By applying adiabatic expansion model we have provided 
a new method to estimate the electron power-law index p 
from single frequency radio observations, which does not 
need correction for the optical depth time delay effects. 
The simultaneous multi-wavelength radio observations of 
GRS 1915+105 reported by Fender et al. (2002) reveal 
short time large variations in the spectral index from op- 
tically thick {a ~ — 0.2)|^ to thin values {a ~ 0.7) which 
corresponds to p of 2.4 in the optically thin regime. They 
observed the fast oscillations of about half an hour period 
on MJD 51751. The power-law index p obtained from our 
observations (Table 1) are in agreement with the measured 
values of Fender et al. (2002) in the optically thin region. 
The simplest explanation to this is that at the peak of 
the oscillation, the contribution from the optically thick 
radio emission is large and hence the delay time correction 
becomes important. However, at low flux level the contri- 
bution from optically thick emission is reduced and the 
delay time correction is less important. This explanation 
for variations in a reported by Fender et al (2002) suggests 
that intrinsically p could be constant which is consistent 
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Fig. 3. Examples of isolated flares from GMRT observa- 
tions. Top panel: The isolated flare seen on 2001 June 28. 
Bottom panel: One of the four minor isolated flares seen 
on July 1. The data is binned at 1 min interval, the time 
is in UT. The dotted line is the fit to the equations (3) 
and (4) for the rising part and for the optically thin decay 
respectively. 

with an adiabatically expanding synchrotron source with 
no additional acceleration. It is difficult to visualise change 
in p over such a short duration of ~ 15 min. 

We can calculate the time delays between emission at 
different wavelengths using this p. The time at which the 
emission at an wavelength A reaches its maximum is given 
by (van der Laan, 1966) 

tm,\ oc (6) 

We calculate the proportionality constant of equation 
6 using our estimate of p and the time delay of 0.25 h 
between 8.3 GHz and 5 GHz radio emission observed in 
GRS 1915+105 by M98. The value of igcm estimated by 
us for May 15, 1997 observations of M98 is ~1.2 h (M98 
gives 0.9 h assuming p — 0) and the ejection of plasma 
occurred at 14.95 — 1.2 = 13.75 h. Inspection of the X- 
ray data at this time (Figure 1 of M98) shows that it still 
falls in the gap in the X-ray data. We have calculated the 
time of formation of the synchrotron bubble for the flares 
observed by us which are given in Table 1. 

We give in Table 2 the observed time delays between 



4 



Ishwarar Chandra, Yadav & Pramesh Rao: Synchrotron Bubbles from GRS 1915+105 



Table 1. Estimated parameters for the synchrotron bub- 
bles from GMRT Observations at 1.28 GHz. 



Obsn 


Time of 


Dur 


Peak 


P 


a 


date 


ejection 


(min) 


Flux 






06/28 


19.8126 


- 30 


38.8 


2.17±0.17 


0.58±0.09 


06/29 


14.0653 


~ 8 


51.3 


1.25±0.12 


0.13±0.06 


06/30 


18.0620 


~ 30 


30.1 


1.64±0.08 


0.32±0.04 


07/01 


15.5077 


~ 35 


27.9 


1.36±0.07 


0.18±0.04 


07/01 


20.3278 


~ 12 


19.4 


1.51±0.13 


0.25±0.06 


07/01 


21.5200 


~ 9 


20.5 


1.12±0.11 


0.06±0.05 


07/01 


21.6818 


~ 6 


32.8 


2.01±0.32 


0.50±0.16 



observations of GRS 1915+105 reported in the literature. 

We have estimated the time delays between these frequen- 
cies using the weighted mean value oi p = 1.46 from our 
observations which are also given in Table 2. For obser- 
vations on M.JD 50675, the calculated time delay of 37 
s for the IR emission is consistent with the observation 
that the spikes in X-ray coincide with the beginning of 
IR flares (Eikcnberry et al. 1998). The time delay of 200- 
400 s reported by them between the X-ray and IR peaks 
is probably due to the time required for the evolution of 
the IR flare as well as for the readjustment of the accre- 
tion disk. Our calculated time delay of 804 s is in good 
agreement with the measured time delay between 15 GHz 
and 8.3 GHz observations reported by M98 (MJD 50583). 
In a detailed study of X-ray - radio observations in GRS 
1915+105 , Klein- Wolt et al. (2002) measure the time de- 
lay of 40 — 45 min from the beginning of state C to the 
peak radio emission at 15 GHz, which agrees with our es- 
timate of 43 min. The observations of M98 (MJD 50700) 
has measured time delay of ^ 960 s between IR and 8.3 
GHz, while our calculated value is 3351 s (0.93h), the dif- 
ference is about 40 min. Inspection of Figure 2 of M98 
show that a previous IR peak is separated by about the 
same amount of time. This suggests that the 8.3 GHz flare 
most likely corresponds to the IR flare ^ 0.9 h prior to 
the 8.3 GHz peak. For observation on MJD 50705, the 
discrepancy of 547 s between the measured and the cal- 
culated time delay between IR and radio emission could 
be due to larger uncertainty in the measured time delay 
because of the presence of two IR peaks. From the radio 
observations at 4.8 and 8.6 GHz on MJD 51751, Fender et 
al. (2002) find a delay time of ^ 600 s based on a broad 
peak (~ 300 s) in the cross correlation function while our 
estimate of time delay is 1030 s. Although the exact an- 
swer to this discrepancy is not clear, it may be due to 
uncertainty in the time delay measurements. 

It may be noted here that the time delays estimated 
using our simple model of synchrotron bubble arc in good 
agreement with many of the observed values between 
different frequencies available in the literature. However, 
naive application of this model predicts much higher 
flux densities than observed at higher frequencies. This 
discrepancy may be the subject of future studies. 



Table 2. Calculated time delays at different frequencies 
for the multiwavelength data from literature. 



MJD 


1^1 


1^2 


M 


At 


Ref 




(Hz) 


(Hz) 


(obs) 


(calc) 




50675 


X-ray 


1.4x10" 


200-400 


37 


E98 


50583 


1.5x10^° 


8.3x10^ 


~750 


804 


M98 


50698 


X-ray 


10^5 


2400-2700* 


2584 


K02 


50700 


1.4x10" 


8.3x10^ 


~960 


3351 


M98 


50705 


1.4x10" 


1.5x10^° 


~2000 


2547 


F98 


51751 


8.6x10^ 


4.8x10'' 


~600 


1030 


F02 



*When calculated from the beginning of State C 
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